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X-ray photoelectron spectra provide a wealth of information on the electronic structure. The extraction of 
molecular details requires adequate theoretical methods, which in case of transition metal complexes has to 
account for effects due to the multi-configurational and spin-mixed nature of the many-electron wave function. 
Here, the Restricted Active Space Self-Consistent Field method including spin-orbit coupling is used to cope 
with this challenge and to calculate valence and core photoelectron spectra. The intensities are estimated 
within the frameworks of the Dyson orbital formalism and the sudden approximation. Thereby, we utilize 
an efficient computational algorithm that is based on a biorthonormal basis transformation. The approach is 
applied to the valence photoionization of the gas phase water molecule and to the core ionization spectrum 
of the [Fe(H 2 0) 6 ] 2+ complex. The results show good agreement with the experimental data obtained in this 
work, whereas the sudden approximation demonstrates distinct deviations from experiments. 

PACS numbers: 31.15.A-, 31.15.aj, 31.15.am, 31.15.vj, 32.30.Rj, 32.80.Aa, 32.80.Fb, 33.60.+q 


I. INTRODUCTION 

Spectroscopy, tracing the changes of the energetic lev¬ 
els upon various physical interactions and in course of 
dynamical processes, is one of the most powerful ana¬ 
lytical tools finding applications in almost every field of 
natural sciences, medicine, and engineering. Among the 
various spectroscopic methods, experiments in the X-ray 
regime are suitable to explore the structure of materi¬ 
als in all aggregation states.® The most popular variants 
of X-ray spectroscopies comprise first order absorption 
and photoelectron spectroscopy as well as second order 
fluorescence and Auger spectroscopy. 1 2 

Upon X-ray irradiation, a core hole is created either 
resonantly or off-resonantly, i.e. an electron is excited to 
a bound or continuum state. Since the binding energies 
of core electrons are significantly different for various el¬ 
ements and core orbitals have a very localized probabil¬ 
ity density, X-ray spectroscopy can be used as an ele¬ 
ment specific local probe of the electronic structure of an 
atom in its environment. This in particular distinguishes 
spectroscopy in the X-ray and UV/VIS regime, since the 
latter gives insight into transitions between delocalized 
molecular orbitals (MO). 3 ^ 

Remarkably, photoelectron spectroscopy was found to 
be sensitive to specific solute-solvent interactions of tran- 
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sition metal (TM) complexes in solutions 4 8 which are 
essential for understanding processes in catalysis, bio¬ 
chemistry, and material sciences PHU However, due to the 
complex electronic structure and notable (especially for 
core ionization) spin-orbit coupling (SOC) of TM com¬ 
pounds, photoelectron spectra (PES) are rich of features, 
and an unambiguous assignment is difficult without the 
aid from theoretical calculations. 1 

On the theory side, a number of methods has been 
proposed for the assignment of PES. The simplest ones 
are Hartree-Fock or Kohn-Sham density functional the- 
Q ry.l 16 l i7 l Here, MO energies are attributed to the PES 
transition energies via Koopmans’ theorem 18 and in¬ 
tensities are not analyzed. In some cases, this may 
provide already a good interpretation of the observed 
spectrum. 5 However, such a simple picture is not able 
to describe more complex effects such as combination 
transitions. For that purpose, methods based on a 
Green’s function approach have been introduced, see e.g. 
Ref. [19j For instance, the algebraic diagrammatic con¬ 
struction formulation 2 ^ enjoys particular popularityP®®^ 
These methods directly deliver the spectroscopic observ¬ 
ables, i.e. transition energies and intensities, as poles and 
residues of the Green’s function, 19 ^ respectively, avoiding 
the calculation of the stationary wave functions. The 
relativistic treatment necessary for TM compounds has 
been implemented within a four-component formalism 24 
which is quite computationally demanding. 

In principle, any quantum chemical method capable of 
describing excited electronic states can be used to obtain 
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the PES peak positions of main and combination transi¬ 
tions as energy differences between the A/'-electron initial 
and N — 1-electron final states. Methods which have 
been previously used in PES calculations include Config¬ 
uration Interaction (Cl) E1H2 9 Time -Dependent Density 
Functional Theory (TDDFT) * 30 * 31 * Equation-of-motion 
Coupled Cluster (EOM-CCSD),®IHl and multiconfigura- 
tional methods based on Complete or Restricted Active 
Space Self-Consistent Field (CASSCF/RASSCF) wave 
function.®® The latter group is of particular impor¬ 
tance for TM compounds, since they are known to some¬ 
times have wave functions which, even in the ground 
state, cannot be represented by a single configuration. 37 
Alternatively, the semi-empirical vale nce bond Ligand- 
Field Multiplet (LFM) technique 1 ® is widely used for 
X-ray photoelectron spectroscopy of TM compounds. 
In all these methods, the SOC is either not included 
or treated within the multi-configurational Dirac-Fock 
method in ^'-coupli ng li mit 39 Jor within semi-empi rically 
parametrized (LFIvtf®) an d a b initio (RASSCF* 35 * 36 h 
LS'-coupling limit. 


II. THEORY 


We consider neutral or charged molecules with N elec¬ 
trons, which are initially in their ground state | In 
the sudden ionization limit HHMS] the final state can be 
written as an antisymmetrized product |^^ _1 '0 el (k)) of 
the continuum state of the photoelectron 2p el (k) and the 
wave function of the N — 1-electron remainder 

Assuming that all photoelectrons with a certain kinetic 
energy (£*. = h 2 k 2 /(2 m G )) are detected regardless of their 
outgoing direction and spin, their number per unit time is 
proportional to the transition rate, including all possible 
initial and final states with energies £j and £p and inte¬ 
grated over all directions, of the outgoing electron. 
In the long wavelength approximation, it reads: 
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The intensities in LFM and most of Cl studies are esti¬ 
mated in the sudden approximation (SA) 40 41 neglecting 
the kinetic energy dependence of the transition strength 
and approximating it in a form of a wave function over¬ 
lap, see Sec. |TT| In principle, the SA should be valid for 
high kinetic energies of the outgoing electron. However, 
for certain applications the SA has been show n to be un¬ 
able to make reliable prediction of intensities! 28 * 34 * 

For TDDFTpnm EOM-CCSD™3] Cl,® and 
CASSCF® based methods the Dyson Orbital (DO) 
formalism^ has been applied. Although computa¬ 
tionally more demanding, it gives a rather accurate 
description of intensities. Further, it provides a compact 
representation of the PES matrix elements by virtue of 
reducing the 7V-particle to a one-particle integration, see 
Sec. |TT] Finally, DOs can be rigorously employed for the 
determination of angular-resolved PES. 43 

The present work sets the focus onto formulating and 
testing a simulation protocol for the incorporation of 
SOC and multi-reference effects into the DO formalism 
for accurate description of PES intensities of TM com¬ 
plexes. Thereby, we employ the RASSCF multi-reference 
approach!®® and include SOC within state interaction 
(RASSI) method 47 in the atomic mean field integral ap¬ 
proximation.® The paper is organized as follows. In 
Sec. [IT] we give a general overview of the theoretical back¬ 
ground and of the developed methodology. Computa¬ 
tional details are provided in Sec. m and experimental 
details in Sec. [IV] Section [V] discusses application of the 
described protocol to two model systems, i.e. gas phase 
water (valence PES) and the [Fe(H20)e] 2+ complex (core 
PES), which corresponds to solvated iron(II) ions in wa¬ 
ter. In Sec. [VT]we conclude that the proposed approach 
allows the description of PES on the same footing as 
processes involving photons, i.e. X-ray absorption and 
resonant inelastic scattering, reported recentlyF 9 56 


where d is the dipole operator, e is the polarization of the 
incoming photon, the lineshape function A(£) accounts 
for the finite width of the excitation pulse, inhomoge¬ 
neous, and other broadening effects. The thermal pop¬ 
ulation of the initial states enters Eq. 0 via the Boltz¬ 
mann factor f B (£i). 

In the following we will specify the many-body states 
fk^ 7 and in terms of Slater determinants (SD), 0^, 

composed of single particle MOs, Let us consider 
the PES matrix elements in Eq. (TTl) for two SDs giving a 
contribution to the initial state , i.e. Oj^, and the fi¬ 
nal state i.e. O^ -1 . Omitting for convenience the 

polarization and k dependence and applying the strong 
orthogonality condition between the free electron func¬ 
tion and initial molecular orbitals = 0? see, 

e.g., Ref. 157) one can write 
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where we introduced the DO 
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Here, V £ Sn denotes all possible permutations of N 
orbital indices with parity p coming from the structure 
of the SDs. Thus, the expression for the PES matrix 
element simplifies to a one-particle integration (electron 
coordinates x&) because the N — 1-dimensional integra¬ 
tion is comprised into the DO with the normalization 
factor VN 
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Note that in general the one-electron orbitals of the initial 
7V-electron and final N — 1-electron bound states are not 
orthonormal due to electronic relaxation upon electron 
removal. The DO can be considered as an analogue of 
the reduced one-electron transition density. Its direct 
correspondence is the wave function of the photoelectron 
before the ionization. Further, the DO is not normalized 
and the PES intensity is proportional to its norm. The 
convenience of the DO formulation of the PES matrix 
elements comes with the fact that the DO contains all 
system- and method-specific information. Therefore, one 
can vary the approximations for ?/+i, without changing 
the DO. 

In practice, the 7V!-fold summation in Eq. © can be 
circumvented by calculating N determinants of TV — 1 x 
N — 1 dimensional overlap matrices, similar to Ref. 0TJ 
Thus the computational effort reduces to 0(N 4 ). A more 
efficient scheme can be obtained upon transformation of 
the nonorthonormal MO basis sets and Cl coefficients of 
the initial and final states and to biorthonor¬ 

mal sets! 47 * 58 ^ This strongly simplifies the DO expression 
in Eq. ©> which reduces to a single term. To show this 
we consider the case where the removed electron was in 
orbital <^J in the initial state If this orbital is not 
contained in the final state O^ -1 , then only one permu¬ 
tation in Eq. |3| leads to nonzero results. That is the 
permutation V, which shifts the electron in (p J into the 
TVth half space of the initial state which does not take 
part in the integration. The parity of the permutation V 
is p = N — r. Thus, we arrive at a simple expression for 
the DO in the biorthonormal basis 

= (-1+“++ (5) 

The most important feature of the representation in 
the biorthonormal basis is that the DO for a pair of de¬ 
terminants corresponds to a spin-orbital of the initial SD 
taken with the appropriate sign. Thus the computational 
effort does not depend on the number of electrons, which 
makes this approach very efficient. 

In the single SD picture, the DO depends on the co¬ 
ordinates of only one electron. Thus, no combination 
transitions, where one electron is ejected and others are 
simultaneously excited, are taken into account. However, 
in a multi-reference description of the DO, these effects 
are included. To obtain wave functions, which account 
for the multi-reference character of transition metal com¬ 
plexes, as well as SOC effects, we follow a two-step strat¬ 
egy. First, the spin-free RASSCF wave functions having 
ground state spin S as well as S' ± 1 are calculated. They 
have the form of a linear combination of different SDs 
©i, weighted with Cl coefficients Ck% (omitting the su¬ 
perscript for the number of electrons in the following) 

Nci 

^ < K=F,I = ( 6 ) 

i= 1 

Assuming that the initial state can be described by a 
single SD (©o), the DO within the Cl approach, omitting 


for convenience / and F indices, can be written in the 
form 
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Here, the brakets denote the N — 1-dimensional inte¬ 
gration, i,j, &,... are indices of the occupied orbitals 
from which the electrons are removed, a, 6,... correspond 
to the unoccupied orbitals into which the electrons are 
moved. The first three sums represent the main lh (sin¬ 
gle hole) and combination 2hlp (two holes, single par¬ 
ticle) and 3h2p (three holes, two particles) transitions, 
respectively. 

The concept of active space within RASSCF as ap¬ 
plied to Eq. ([t]) allows to flexibly vary the highest order 
of correlation terms, in principle up to full Cl within the 
given orbital subspace. However, this subspace needs to 
be large enough to include all relevant ionization chan¬ 
nels. Interestingly, when employing the biorthonormal 
basis set transformation, orbital relaxation effects, which 
can cause satellites due to non-orthogonality of initial 
and final orbitald^, are shifted completely to correlation 
effects, i.e. the Cl expansion. 

At the second step, these wave functions are coupled 
with the SOC operator in atomic mean-field integral 
approximation 48 to generat e spin -orbit wave functions 
within the RASSI approach. 47 5 ^* These spin-orbit wave 
functions are expanded in terms of the spin-free states 
wh ere K = F, /, (Eq. 0) with spin cr ( N a 
states in total), magnetic quantum numbers M a and SOC 
coefficients 


5+1 N a a 

^W=F,I = ^ 

cr=S— 1 n= 1 M a =—<j 


(8) 


The particular choice of spin manifolds is dictated by the 
SOC selection rules AS = 0, ±1. 

Thus, the multi-configurational DO including SOC for 
the transition from initial state I to final state F can 
be expanded in terms of single determinant DO T 8D 
(Eq. §) 
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Here, the nested sums over spin, magnetic quantum num¬ 
ber and spin-free states are replaced by one sum over all 
different spin-orbit states with the total number of final 
and initial SOC states TVfP and Nf°, respectively. The 
above Eq. © can be used to calculate the PES matrix 
elements = ^'0 eZ |d|<I>|O^, which is the main working 
expression used in this work. 

Since an electron with either a or /3 spin is removed 
from the initial spin-free states with spins S and S± 1, the 
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Figure 1. Schematic of the directly spin-coupled manifolds 
of states (the double-sided arrows, AS = 0, ±1) of the N and 
N— 1 electron systems. Depending on the spin of the outgoing 
electron, the spin of final N — 1 electron states changes by 
AS = ± 1/2 (single-sided arrows). 


final states with AS = ±1/2 need to be considered. That 
is why in general seven spin manifolds of the unionized 
molecule with N electrons and its N — 1-electron ion need 
to be taken into account as depicted in Fig. [l] 

Frequently, the SA is used for the calculation of the 
transition dipole matrix element in Eq. Q SSED Here, 
the matrix element is approximated as an overlap integral 
neglecting the k-dependence of the transition strength 
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where the operator a annihilates the electron from the 
occupied MO. Expressed in terms of DOs in the SA ap¬ 
proach, the one-electron integration for (^ el |d|$.pj) is 
omitted and the intensity is approximated by the DO 
norm, |<Ff/| 2 , only. Most of the published papers use the 
SA for the prediction of the intensities of the combina¬ 
tion transitions relative to the intensities of main lines. 
However, the accuracy of the predicted intensity ratios 
between different main transitions occuring at high exci¬ 
tation energies has been questioned in Refs. 1281591 


III. COMPUTATIONAL DETAILS 

The spectra were calculated for two test systems, i.e. 
gas phase water and the [Fe(H 20 )e] 2+ complex mimick¬ 
ing the Fe 2+ ion within its first solvation shell in aqueous 
solution. Thus the following processes were studied 

H 2 0° —> H 2 0+ + e- 

[Fe(H 2 0) 6 ] 2+ —► [Fe(H 2 0) 6 ] 3+ + e" . 

The geometry of water was optimized at the density 
functional theory level with the B3LYP functional 60 61 


together with the 6-311G(d) basis set! 62 * 63 1 For the 
[Fe(H 20 ) 6 ] 2+ complex the geometry was first opti¬ 
mized using the MP2 method with the cc-pVTZ ba¬ 
sis set!®^ Then the Fe — O bond lengths were set to 
the CASPT2(10e _ , 12MO)/ANO-RCC values obtained 
in Ref. [671 All geometry optimizations were done with 
Gaussian 09 program packaged 

For state-averaged RASSCF calculations of water, the 
full valence RAS2 active space additionally including the 
Is core orbital of oxygen was chosen (Fig. [2]). For the 
iron complex, the active space included the 3d orbitals of 
iron in the RAS2 space and allowed single electron exci¬ 
tations from the iron 2p core orbitals in the RAS1 space. 
This results in the singly core-excited or ionized states, 
having full flexibility in the valence 3d manifold. This 
active space corresponds to that used for the study of X- 
ray absorption and resonant inelastic scattering spectra— 
providing interpretation on the same footing. The wave 
function was first optimized for the ground state and then 
all orbitals apart from the active ones were kept frozen. 

The (8s4p3d)/[3s2pld] ANO-RCC basis on hydro¬ 
gen and (14s9p4d3f)/[5s4p3d2f] basis on oxygen— 
were employed for water, which corresponds to 
quadruple-zeta quality. For Fe 2+ the contractions 
(21sl5pl0d6f)/[6s5p3d2f], (14s9p4d3f)/[4s3p2dlf], and 
(8s4p)/[2slp] were used for iron, oxygen and hydrogen, 
respectively, corresponding to the triple-zeta level. 

To account for dynamic correlation, the second 
order perturbation theory correction (RASPT2) was 
calculated 70 71 for the case of water. To avoid intruder 
state singularities the imaginary level shift of 0.4 Hartree 
was applied. Because of the much larger number of tran¬ 
sitions only RASSCF energies were calculated for the iron 
complex. 

To incorporate scalar r elativ istic effects, the Douglas- 
Kroll-Hess transformation 72 73 up to second order was 
utilized. For the case of water, SOC was not taken 
into account and only transitions from the singlet ground 
state to the 490 doublet singly ionized final states were 
considered. This corresponds to one out of six possible 
branches in Fig. [l] In case of the iron complex, quintet 
(S = 2) and triplet (5 = 1) initial states as well as sextet 
(5 = 5/2), quartet (5 = 3/2), and doublet (5 = 1/2) 
final states were taken into account. Septet (5 = 3) 
and octet (5 = 7/2) states are not possible with the ac¬ 
tive space chosen here. However, it was shown that they 
play a very minor role for X-ray absorption spectrum^ 
and can be neglected. Thus, only five branches out 
of the seven in Fig. [l] were considered. In total, for 
[Fe(H 20 )e] 2+ one initial ground and 1260 core-excited 
SOC states were included. Thermal population of the 
low lying initial states was neglected. 

The RASSCF/RASPT2/RASSI calculations were 
done without any symmetry restriction using a locally 
modified MOLCAS 8.0 program package.^ 

The transition dipole matrix elements in Eq. © were 
calculated with the ezDyson 3.0 program 75 ’via numerical 
integration of the DO with the free electron wave func- 
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Figure 2. Active spaces used for RASSCF calculations of (a) water and (b) [Fe(H20)e] 2+ - 


tion ^(k). Neglecting the interaction between the pho¬ 
toelectron and ionic remainder one may express ^ el (k) 
as a plane wave expanded in a basis of spherical waves 
with spherical harmonics Y/ ?m (k) in k-space as coeffi¬ 
cients (see, e.g., Ref. EH- Further, the spherical waves are 
expanded in position space using spherical Bessel func¬ 
tions ji(k • r), where k = |k|, r = |r|, and spherical har- 
monies Y Lm (r) yielding 

OO l I ^ 

Y(k) =EE • r)Y ltm (r)Yi* m (k) . (11) 

z—✓ z—/ Y 7J- 

Z=0 m=—l 


This particular form of the plane wave has the advantage 
that one can truncate the infinite angular momentum ex¬ 
pansion at some / max , e.g., according to the dipole selec¬ 
tion rule Al = ±1. The natural choice of maximum an¬ 
gular momentum in the expansion Eq. (11) corresponds 
to the / max = Zbas + 1, where /bas is maximum angular 
momentum included in the atomic basis set. However, 
due to lower than spherical symmetry of the molecules 
under study, sometimes a larger / max needs to be selected, 
see Ref. l32land Sect. [V| In case of water, / max = 7, box 
with a side length of 10 A and an equidistant grid of 
300 x 300 x 300 points for numerical integration of Eq. © 
ensures convergence of the DO norms with an accuracy 
of 10 -5 . Since the intensity of transition scales quadrati- 
cally with the norm of the DO, they have been evaluated 
only, if the norm of the respective DO was larger than 
10 -3 . For the iron ion, only contributions from those 
DOs that have a norm larger than 10 -2 were taken into 
account. Here Z max = 5, box size of 16A and a grid of 
480 x 480 x 480 points was used for the numerical inte¬ 
gration reproducing the norms of DOs with an accuracy 


of 10- 4 . 

As shown in Sec. [TTJ the spin-coupled DO comprises 
contributions from different spin-states. Therefore, the 
total spin of the DO and outgoing electron is not well 
defined. The consequence is that the DO consists of both 
a and (3 spin wave functions 

$19 = *??(<*) +$19(0), (i2) 


which are complex-valued due to the SOC. Here, we ne¬ 
glect spin-coherence and calculate the squared PES ma¬ 
trix element as 
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where Re and Im represent the real and imaginary parts 
of Dpi. The broadenings of the spectral lines were fit¬ 
ted to reproduce the experimental data and are discussed 
in the next section. No nuclear vibrational effects were 
taken into account. The PES matrix element was inte¬ 
grated over all possible outgoing directions of the pho¬ 
toelectron and averaged over all possible orientations of 
the molecules to mimic free tumbling of the solute in a 
liquid phase. 


IV. EXPERIMENTAL DETAILS 

The PES were measured from a 15 /im vacuum liquid 
jet ! 77 ! 78 1 at the soft-X-ray U41 PGM undulator beamline 
of the Berlin synchrotron radiation facility, BESSY II, 
Berlin. The jet velocity was approximately 100 m • s _1 , 
and the jet temperature was 6°C. At operation con¬ 
ditions, the pressure in the interaction chamber is ~ 
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1.5 • 10 -4 mbar. Electrons were detected normal to both 
the synchrotron-light polarization vector and the flow of 
the liquid jet. A 100 fi m diameter orifice that forms the 
entrance to the hemispherical electron energy-analyzer 
(Specs Leybold EA10) is typically at approximately 0.5 
mm distance from the liquid jet. Because of the small 
focal size (12 x 23 fim 2 ) of the incident photon beam at 
the interaction point with the liquid jet, by moving the 
jet slightly out of focus the photoelectron spectrum from 
gas-phase water can be measured. The valence spec¬ 
trum of the latter was obtained using 180 eV photon 
energy; the total energy resolution was better than 60 
meV. For the aqueous phase Fe 2+ 2p core-level spectrum 
the photon energy was 925 eV, and the energy resolu¬ 
tion was approximately 300 meV. The 2M iron aqueous 
solution (pH « 2.5) was prepared by adding anhydrous 
FeCl 2 salt (98% purity, Sigma Aldrich) to highly dem¬ 
ineralized water (> 17 MOhm • cm -1 ) water. The ma¬ 
jor species in this pale-lime colored aqueous solution are 
[Fe(H 2 0) 5 Cl]+ (« 77%) and [Fe(H 2 0) 6 ] 2+ (« 23%)P 


V. RESULTS AND DISCUSSION 
A. Water 

The PES of water was chosen as a benchmark for the 
derived protocol for calculation of PES for two reasons. 
First, the water valence PES is well und ersto od both 
experimentally and theoretically ! 25 ! 26 ! 28 * 59 * 80 *^ Second, 
water is not influenced by strong SOC and has mostly 
single-configurational character as long as only the low¬ 
est valence excitations are regarded. This substantially 
simplifies the calculation and analysis of the nature of 
transitions, since DOs represent mostly single MOs. 

In Fig. [3j the experimental and theoretical PES for 
the 180 eV energy of the incoming photon are shown. 
For convenience, the individual transitions are also given 
as a stick spectrum. The normalized intensity is plotted 
against binding energy E b = £p — £j. where £p and £j 
are the energies of final and initial states, respectively. 
Additionally, the DOs corresponding to selected transi¬ 
tions are presented. The broadening of each transition 
was fitted to reproduce the experiments: peak (1) (see 
Fig. [3| Lorentzian lineshape FWHM=0.17 eV, peaks (2) 
and (3) Gaussian lineshape with FWHM of 1.18 eV and 
1.75 eV, respectively. For all other transitions contribut¬ 
ing to peak (4), Gaussian profiles with widths 1.5 eV for 
E b < 35 eV and 3.0 eV for E b > 35 eV were applied. The 
absolute energy shift of -0.78 eV was chosen to align peak 
(1) with the experiment. 

The calculated spectrum is in rather good agreement 
with the experimental data. There are slight variations 
of intensities of peaks (2) and (3) relative to (1) and the 
most notable discrepancies between theory and experi¬ 
ment are observed for peak (4). These deviations in the 
intensities and lineshapes could be ascribed to the fact 
that nuclear vibrational effects were not taken into ac- 
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Figure 3. Calculated and experimental PES of water in 
the gas phase for 180 eV excitation energy. Full calculation 
corresponds to numerical integration of PES matrix element, 
SA means sudden approximation.^^ Spin-free DOs of the 
selected transitions are also shown. 


count. However, the relative energetic positions of the 
peaks are predicted with very high accuracy. To note 
is the fact that RASPT2 correction is essential here to 
reproduce the transition energies. This stems from the 
fact that in RASSCF dynamic correlation is accounted 
for in an unbalanced way within the active space and 
may substantially change upon removal of one electron. 
To correct for this behavior a more complete treatment 
of correlation, such as RASPT2, is needed. 

It can be seen that SA predicts the relative intensities 
of peaks (1) to (3) with similar quality as the full cal¬ 
culation employing integration of the DO with the free 
electron wave function. But, for peak (4) the agreement 
between SA and experiment is not even qualitative. 

The assignme nt of peaks fully agrees with that estab¬ 
lished previously!^* 26 f 28 * 59 * 8Q * 811 Peaks (l)-(3) correspond 
to ionizations from the lone pairs no(l&i) and no(3ai) 
and the cr p (l& 2 ) MO of water, respectively. The fourth 
peak consists of several transitions of different nature 
most of them having <j s (2ai) character. This can be as¬ 
cribed to the multi-configurational character of the wave 
function and appearance of combination transitions, see 
Supporting Information. Interestingly, the DOs of a\ 
symmetry resemble very closely the ground state Hartree- 
Fock MOs, but in fact they are linear combinations of the 
cr p (2ai) and cr p (3ai) RASSCF active orbitals (cf. Fig.[2|. 

The water molecule has rather low point symmetry 
and convergence of the series in Eq. ( [IT] ) is quite slow. 
The contributions of the partial {/,m|waves to inten¬ 
sities of peaks (1) and (3) are presented in Fig. [4j The 
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DOs of transitions ( 1 ) and ( 2 ) represent almost pure 2 p 
orbitals of oxygen and in accordance with the dipole se¬ 
lection rules (A l = ± 1 ) only contributions with l = { 0 , 2 } 
notably differ from zero. In contrast, the DOs for tran¬ 
sitions (3) and (4) deviate strongly from atomic charac¬ 
ter, being rather delocalized MOs such that the series 
in Eq. © converges notably slower, approaching zero 
only for Z > 6 . In general, one can conclude that the 
more the DO is delocalized over the molecule and the 
less symmetric it is, the more terms have to be included 
in the photoelectron wave expansion (Eq. ([II))), see also 
Ref. M 


B. Fe 2+ (aq.) 

The calculated L-edge core PES of [Fe(H 2 0)6] 2+ at 
incoming photon energy of 925 eV is shown in Fig. 
together with experimental results for a 2M aqueous so¬ 
lution of FeCl 2 at the same energy. The stick spectrum 
was broadened using the Voigt profile 

/ oo 

G(cr, x')L( 7 , x — x')dx , (14) 

-oo 

with the Gaussian and Lorentzian lineshape functions 
G(cr, x) and L(y, x). For the broadening of the L 3 peak 
(£5 > —727.4 eV), 0.5 eV and 0.7 eV were used for 
the Lorentzian and Gaussian FWHM, respectively. For 
the I /2 peak (below E^ < —727.4 eV), 0.7 eV for both 
Lorentzian and Gaussian width in the Voigt profile was 
used. Additionally, the calculated spectrum was shifted 
as a whole by +10.65 eV for better comparison with the 
experimental data. 

The spectrum consists of two prominent bands called 
L 3 and +2 which correspond to the SOC components of 
the created core-hole, i.e. 3/2 and 1/2 total angular mo¬ 
mentum, respectively. Transitions to states 253 — 1000 
and 1001 — 1260 of the [Fe(H 20 ) 6 ] 3+ ion correspond to 
the L 3 and L 2 peaks, respectively, in the core PES shown 
in Fig. [5] The computed core PES with numerical inte¬ 
gration according to Eq. © (full calculation) is in rather 
good agreement with the experiment. Our method re¬ 
produces well the T 2 /I /3 energy splitting, intensity ra¬ 
tio and the general asymmetric shape of the bands with 
the long tails at the low energy sides. Almost all of the 
1260 transitions have notable intensity, hence the bands 
and corresponding tails are formed by hundreds of lines, 
most of them being combination transitions. The small 
discrepancies, e.g., the small peak at —722.9 eV and the 
minimum at —720.5 eV, also originate from a large num¬ 
ber of transitions of different origin, what hinders a de¬ 
tailed analysis. In principle, they could be ascribed to 
the lack of dynamic correlation (no RASPT 2 correction) 
and to the presence of different species in solution. Note 
that ionization only from [Fe(H 2 0)6] 2+ was considered in 
computations. Further, the averaging over several ther¬ 
mally populated electronic initial states could be neces¬ 
sary. Finally, due to the high density of states, the tails 



Figure 4. Normalized contributions of different {/, m} partial 
waves to the intensity of selected transitions: (a) and (b) tran¬ 
sitions (D and (3) of an isolated water molecule, respectively, 
see Fig.Dj (c) real part of ft spin DO contribution of transition 
(2) of the [Fe(H20)e] 2+ core PES, see Fig. [H] Each group of 
21 + 1 sticks corresponds to m- components, see Eq. (11). Cor¬ 
responding DOs are shown as insets using different contour 
values for visual clarity. 


of the +2 and £3 peaks are very sensitive to the wings of 
the lineshape function. Therefore, an inclusion of more 
than two different sets of broadening parameters might 
be necessary as well. 

The wave functions of the final states of [Fe(H 20 ) 6 ] 3+ 
are much more complex in comparison to those of wa¬ 
ter. First of all, the core-ionized final states do not have 
a leading contribution from one configuration and repre¬ 
sent a combination of many configurations with compara¬ 
ble weights in Eq. ©• To illustrate this issue, the weights 
of the three most important configuration state functions 

























B Fe 2+ (aq.) 



( 1 ) ( 2 ) ( 3 ) 



Figure 5. (a): Experimental (2M FeCE aqueous solution) 

and calculated (for [Fe(H20)e] 2+ cluster) core PES for incom¬ 
ing photon energy of 925 eV. Full calculation corresponds to 
numerical integrati on of PES matrix element, SA means sud¬ 
den approximation. 40 41 (b): Real and imaginary parts of a 
and (3 spin contributions to the DOs for selected transitions. 


for the quartet final states of [Fe(H 2 0)e] 3+ are shown in 
the SI. For other spin manifolds the dependences look 
similar and are not shown. Second, due to strong SOC 
for the core-ionized states, a pronounced spin-mixing of 
states within the LS'-coupling scheme applied here is ob¬ 
served. More details of this mixing are given in the SI. 
This implies that the spin of the final states is not well 
defined and their wave functions represent linear com¬ 
binations of the sextet, quartet, and doublet spin-states 
(Eq. ©). The complex structure of the wave function 
is reflected in the DOs (Fig. [5 ]d) . Here, most transitions 
correspond to DOs, where the real and imaginary parts 
of the a and /3 spin contributions are of comparable mag¬ 
nitude. Thus a calculation of only real DOs with definite 
spin would produce erroneous intensities. An exception 
occurs at the rising flank of L 3 (DO (3) in Fig.^) which 
can be ascribed to nearly pure quintet-sextet transitions 
and thus corresponds to photoelectrons with f3 spin. The 
/3 electron is easily removed because of the pairing (ex¬ 


change) energy. 

The DOs for different transitions of both the and 
L 3 bands are combinations of pure atomic 2p orbitals of 
iron and differ in the absolute and relative magnitude 
of real and imaginary a and f3 components. This corre¬ 
sponds to the fact that the electron is removed from the 
core 2p orbitals and, in contrast to water, the series in 
Eq. © converge to zero very quickly for all transitions. 
One example is shown in Fig. |§ for the most intense L 3 
transition. Since only atomic 2pp e orbitals are contribut¬ 
ing to the DOs, the dipole selection rules hold strictly. 

In Fig. [5] we also show the SA results, where the inten¬ 
sities of the transitions are approximated as the norms of 
DO and thus the computationally demanding integration 
in Eq. © is avoided. Apparently, the SA gives a PES 
which substantially deviates from both the experiment 
and the full calculation, most notably it shows rather in¬ 
tense tails of Z/ 3 and L 2 bands. This deviations cannot 
be eliminated by the fitting of broadenings; note that in 
Fig.[5]the same broadening parameters as for the full cal¬ 
culation are used. This result shows th at alt hough this 
approximation is intensively used in Cp® 9 and LFM 
calculations of core PE^and in RASSCCF/RASSI cal- 
culations of valence PES of heavy elements p- 6 it cannot 
be considered as being generally accurate. This holds 
especially if the kinetic energy of the ejected electron is 
relatively low. Thus to obtain accurate L-edge core PES 
the integration of the DO with the free electron wave 
function is required. Commonly, the SA is used when 
the photoelectron momentum angular distribution is not 
of interest. Our results point to the failure of the SA 
even in such cases. 


VI. CONCLUSIONS 

We have presented a multi-reference approach to core 
and valence photoelectron spectra of transition metal 
complexes, taking into account the essential effects due 
to the multi-configurational character of the wave func¬ 
tion and spin-orbit coupling. This method is an exten¬ 
sion of the D vson orbital formalism, previously applied 
with TDDFT™ an d EOM-CCSeSSEH techniques, to 
RASSCF/RASSI wave functions. Thereby, an essential 
point for efficient computation of the Dyson orbitals is 
the biorthonormal MO basis transformation. 58 The pro¬ 
posed protocol includes the numerical integration of the 
matrix elements, resulting in more reliable intensities as 
compared with the widely used sudden approximation. 

This present approach is complementary to the theo¬ 
retical X-ray spectroscopic techniques, which have been 
recently reported for absorption, fluorescence, and inelas¬ 
tic scattering spectraP^IO Having the method describing 
photon-out and electron-out events on the same footing 
provides additional tools to address the electronic struc¬ 
ture of transition metal complexes. 

The computational protocol has been demonstrated 
for two examples, i.e. the valence PES of gaseous wa- 
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ter and L-edge core PES of the [Fe(H 2 0)6] 2+ cluster 
as a model for the Fe 2+ ion in aqueous solution. In 
both cases the agreement between theory and experi¬ 
ment was rather good. In particular for the aqueous ion 
the RASSCF/RASSI wave function has a rather com¬ 
plex structure. This is reflected in the DO, i.e. due 
to spin-mixing the Dyson orbitals for core-ionization of 
[Fe(H 2 0)e] 2+ contain complex-valued a and /? contribu¬ 
tions. This immediately implies that the usage of real 
DOs with definite spin would give erroneous spectra. 

We have contrasted our results with those of the widely 
used sudden approximation. The latter shows for rela¬ 
tively low excitation energies notable deviations from ex¬ 
perimental spectra for the aqueous Fe 2+ ion. Hence, in 
general the numerical integration of the Dyson orbitals 
with the free electron wave function should not be dis¬ 
missed. 
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